Superconducting photonic crystals 
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The band structure of a novel type of photonic crystal with superconducting constituent elements 
is calculated numerically via a plane wave expansion. The density of states and the dependence of 
the width of the photonic gap on the filling factor is analyzed for a two-dimensional photonic crystal 
consisting of an infinite array of parallel superconducting cylinders. 
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Photonic crystals, artificial media with a spatially pe- 
riodical dielectric function that were first discussed by 
Yablonovitch— and John£, are the subjects of growing 
interest due to various modern applications^^ This pe- 
riodicity can be achieved by embedding a periodic ar- 
ray of constituent elements ("particles") with dielectric 
constant e\ in a background medium characterized by 
dielectric constant 62- Different materials have been 
used for these elements: dielectrics, semiconductors and 
metals, 5 ! 6 ! 7 ! 8 ! 9 ! 10 

Photonic gaps are formed at frequencies u> for suffi- 
ciently high dielectric contrast w 2 (ei{u) — £2(0;)), be- 
cause the quantity vj 2 e(u)) enters in the wave equation. 5,6 
Thus, only metal-containing photonic crystals can main- 
tain the necessary dielectric contrast at small frequencies 
due to Drude- like behavior EMet(w) ~ — l/u; 2 £iilBut the 
damping of electromagnetic waves in metals can suppress 
many useful properties of metallic photonic crystals. 

We suggest in this Letter a novel type of photonic 
crystal consisting of superconducting elements (embed- 
ded in dielectric media) or inverse superconducting crys- 
tal with (vacuum) holes. This photonic crystal can be 
used for (THz) frequencies smaller than superconducting 
gap where the damping is negligible. 

For simplicity, we consider superconducting particles in 
the London approximation (i.e. for London penetration 
length Sl = [mc 2 /(Aime 2 )] 1 / 2 £, £ is the coherence 
length). The optical properties of superconductors have 
been extensively studied (see, for example, Refs. fllT^ 
and references cited therein). At (T c — T)/T c <C 1 and 
hut <C A <C T c a simple relation relation for the current 
density is given bjiii 
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where cr is a conductivity of normal metal. 
Substituting Eq. Q into the wave equation 
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and seeking solutions with harmonic time variation of 
the electric field, i.e., E(r,t) = E (r)e lwt , E = iw/cA, 



we have for a two-component photonic crystal: 
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where rj(r £ S) = 1 if r is inside superconductor elements 
S, and otherwise rj(r £ S) = 0, e s (oj) being the dielectric 
function of the superconductor and Eq the dielectric con- 
stant of the surrounding medium. Here we adopt the ap- 
proximation e s (uj) — l+c 2 /(<5 2 oj 2 ), which neglects damp- 
ing in the superconductor. The summation in Eq. PJl 

goes over all lattice nodes a{nj} = Y2i=i n i a i charac- 
terizing positions of superconducting particles (a{n.;} is 
the coordinate space vector of a superconducting particle 
occupying a particular node of the lattice) . 

Here we calculate the band structure of a two- 
dimensional photonic crystal consisting of an infinite ar- 
ray of superconducting cylinders that form a square lat- 
tice. Specifically, by expanding the electric field in plane 
waves (the plane wave expansion method commonly used 
for band structure calculations of solids^ was applied for 
calculation of the band structure of 3D photonic crystal 
with dielectric spheres in the nodes in Ref. [1 4| ) we reduce 
the Maxwell Equation (3) for a superconducting photonic 
crystal to the following eigenvalue problem: 



( <W(k + G) 2 - ^M GG , ) E k (G') 
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where P = c/5l\ G = iiGi + 12G2 is the reciprocal 
lattice vector for 2D square lattice. Here i\ and %2 are 
arbitrary integers, and Gi and G2 are 2D orthogonal vec- 
tors: if the unit cell vectors ai and a2 of the 2D square 
lattice are oriented along x and y real space unit vec- 
tors, then Gi = 2nx/a and G2 = 2iry/a, with a being 
the lattice spacing. Furthermore, Ek(G) are the Fourier 
components of the electric field eigenfunction E k (x) . The 
matrix Mgg' is related to the Fourier transform of the 
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periodic dielectric function in this superconducting pho- 
tonic crystal: 

e(G-G') = 5 GG , + ^M GG ,, (5) 

UJ Z 

where 

M GG , = fe + (f-l), G = G' 

M GG > = 2/(e-l) Jl | ( G G _' G y ) , G^G'. (6) 

Here / = S supercon d/ S = irr 2 /a 2 is the filling factor of 
the superconductor(S , supercon( i is the cross-sectional area 
of the cylinder [in the the plane perpendicular to the 
cylinder axis]; S is the the total area occupied by the 
lattice unit cell, r is the cylinder radius, a is the lattice 
spacing); finally, J\ is a Bessel function. 

The band structure for the filling factor / = 0.3 
and P = 5 x 2nc/a is presented on Fig.l (we use as 
the unit of frequency the lattice frequency 2irc/a)}^- 
For a YBCO superconductor 6l ~ 200nm, and thus 
P 1.499 x 10 15 s -1 . We have picked the value of P = 5 
in units of 2nc/a to attain a significant band gap in the 
spectrum. Since for YBCO the value of P is fixed, our 
results with P — 5 x 2irc/a correspond to the specific 
lattice spacing a = 6.283/im. The density of photonic 
states for the crystal presented in Fig. 2 is defined as the 
total number of different wavevectors k in all directions 
in one unit cell of the reciprocal lattice corresponding to 
the dispersion curve uj = ui{k) in a small frequency region 
divided by the frequency range spanned. As illustrated 
in Fig. 3, the width of the photonic gap is a nonmonotonic 



function of the filling factor /. A photonic gap exists only 
in the range 0.23 < / < 0.44. The widest gap (A = 0.35) 
corresponds to / = 0.32. 

In summary, we have obtained the band structure of a 
superconducting photonic crystal. The width of the pho- 
tonic gap for the system considered is a non-monotonic 
function of the filling factor. The advantage of a pho- 
tonic crystal with superconducting particles is that the 
dissipation of the incident electromagnetic wave due to 
the imaginary part of the dielectric function is much 
greater for normal metallic than for superconducting par- 
ticles, because the imaginary part of the dielectric func- 
tion for superconducting particles is negligible in compar- 
ison with the imaginary part of the dielectric function for 
normal metal particles at frequencies smaller than the su- 
perconducting gap. Thus, in this frequency regime, for 
a photonic crystal consisting of several layers of metallic 
scatterers the dissipation of the incident electromagnetic 
wave by an array of superconducting particles is expected 
to be less than that obtained from an analogous array 
composed of normal metallic particles. Note that exper- 
imental studies of superconducting metamaterials have 
been carried out recently^ 
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FIG. 1: Dispersion dependence for a two-dimensional super- 
conducting photonic crystal with square lattice consisting of 
infinite cylinders having circular cross section. E-polarization 
is considered, with / = 0.3 and P — 5 x 2irc/a. The ordinate 
plots frequencies in lattice units 2nc/a. A band gap is clearly 
apparent. 
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FIG. 2: The density of states g(u>) for a two-dimensional 
superconducting photonic crystal with square lattice (E- 
polarization) with the filling factor / = 0.3 and P = 5x2irc/a. 
The abscissa plots frequencies in lattice units 2%c/a. 



■t— ' 

03 
13 




FIG. 3: Dependence of width of the photonic band on the 
filling factor for a two-dimensional superconducting photonic 
crystal with square lattice in case of E-polarization (P = 5 x 
2iTc/a). The ordinate plots frequencies in lattice units 2nc/a. 



